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ABSTRACT 

Sk 183 is the visually-brightest star in the N90 nebula, a young star-forming region in the Wing of the 
Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope 
which reveals Sk 183 to be one of the most massive O-type stars in the SMC. Classified as an 03-type 
dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption which 
suggests a later type. We propose that Sk 183 has a composite spectrum and that it is similar to 
another star in the SMC, MPG 324. This brings the number of rare 02- and 03-type stars known in 
the whole of the SMC to a mere four. We estimate physical parameters for Sk 183 from analysis of its 
spectrum. For a single-star model, we estimate an effective temperature of 46±2kK, a low mass-loss 
rate of ^ 10~ 7 M yr -1 , and a spectroscopic mass of 46^g Mq (for an adopted distance modulus of 
18.7 mag to the young population in the SMC Wing). An illustrative binary model requires a slightly 
hotter temperature (~47.5kK) for the primary component. In either scenario, Skl83 is the earliest- 
type star known in N90 and will therefore be the dominant source of hydrogen-ionising photons. This 
suggests Sk 183 is the primary influence on the star formation along the inner edge of the nebula. 
Subject headings: open clusters and associations: individual (N90; NGC 602) — stars: early-type - 
stars: fundamental parameters — stars: individual (Sanduleak 183) 



1. INTRODUCTION 

The Magellanic Clouds provide an excellent example 
of interacting, metal-poor dwarf galaxies on our own 
doorstep, with neutral hydrogen maps revealing large 
structures such as the Magellanic Bridge (which connects 
the Clouds), and the extensive Magellanic Stream and 
Leading Arm (e.g., Putman 2000). The Bridge is linked 
to the eastern side of the Small Magellanic Cloud (SMC) 
by an extended 'Wing', which has a lower content of gas, 
dust and stars than the main body of the SMC, suggest- 
ing it should be unfavourable to star formation. How- 
ever, the young stellar cluster NGC 602 presents spec- 
tacular evidence to the contrary, with a star-formation 
rate which appears comparable to well-studied Galactic 
regions (e.g. Cignoni et al. 2009). 

NGC 602 is a loose term for a rich region compris- 
ing at least three distinct clusters (Westerlund, 1964). 
NGC 60243 is embedded in the N90 emission nebula 
(Hcnizc, 1956), NGC 602b is a closely-related, smaller 
cluster approximately l'.5 to the north, and NGC 602c is 
~11' to the northeast. Other nearby Ha emission regions 
include N88 and N89 (Lindsay 104) which, together with 
NGC 602, are associated with the only Ha supergiant 
shell known in the SMC ('SGS-SMC1', Meaburn, 1980). 

1 UK Astronomy Technology Centre, Royal Observatory Ed- 
inburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

2 Institute for Physics and Astronomy, University of Potsdam, 
14476 Potsdam, Germany 

3 Department of Astronomy, University of Wisconsin- 
Madison, 5534 Sterling, 475 North Charter St., Madison, WI 
53706, USA 

4 Astronomy Department, University of Illinois, 1002 W. 
Green Street, Urbana, IL 61801, USA 

5 Scottish Universities Physics Alliance (SUPA), Institute for 
Astronomy, University of Edinburgh, Royal Observatory Edin- 
burgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

6 Cluster 105 from Lindsay (1958) 



Optical imaging of NGC 602a from the Hubble Space 
Telescope (HST) Advanced Camera for Surveys (ACS) is 
particularly striking 7 ] engendering the popular moniker 
of the 'Hubble Bubble'. The HST observations and deep 
infrared images from the Spitzer Space Telescope have 
given us an exquisite view of the stellar content and star 
formation in NGC 602a and N90 (Carlson et al. 2007, 
2011; Gouliermis et al. 2007; Schmalzl et al. 2008; Nigra 
et al. 2008; Cignoni et al. 2009). 

As part of a new study of SGS-SMC1 we have obtained 
optical spectroscopy of its luminous stellar population 
with the Very Large Telescope (VLT). In the course of 
the program we observed Sk 183 (Sanduleak, 1969^ the 
visually brightest star in the N90 nebula (see Figure [l]). 
In this article we present the classification and analysis 
of these observations, which reveal it to be one of the 
most massive O-type stars in the SMC. 

2. SPECTROSCOPY 

Optical spectroscopy was obtained on 2010 October 
24-26 with the VLT FLAMES instrument (Pasquini et 
al. 2002), as summarised in Table [T] Skl83 was ob- 
served as part of the larger survey of the region with the 
Medusa- fibre mode of FLAMES, using three of the stan- 
dard settings of the Giraffe spectrograph: LR02, LR03, 
and HR15N, delivering spectral resolving powers of 7000, 
8500, and 16 000, respectively (e.g. Evans et al. 2011). 
Individual exposure times were 1800s for all settings. 

The ESO Common Pipeline Library FLAMES reduc- 
tion routines (v. 2. 8. 7) were used for the standard process- 
ing stages, i.e. bias subtraction, fibre location, summed 
extractions, division by a normalized flat-field, and wave- 

7 http: //hubblesite.org/ newscenter/archive / releases/2007/ 04 / 

8 a = 01 h 29 m 24! 6, <5 = -73° 33' 16'.'43, J2000. Aliases include: 
NGC 602a #15 (Westerlund, 1964); NGC 602 #8 (Hutchings et al. 
1991); SMC 83235 (Massey, 2002). 
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Table 1 

Summary of VLT-FLAMES spectroscopy of Sk 183 



Figure 1. HST-ACS F814W image of NGC602a and its asso- 
ciated H II region (N90), centred on Skl83. NGC602b is ap- 
proximately 1'5 to the north. The 10 pc s cale assumes the 'short' 
distance modulus of 18.7 mag (see Section |3.1[ l. 

length calibration. The spectra were then corrected to 
the heliocentric frame, and a median sky spectrum was 
subtracted (for more details see Evans et al. 2011). 

Qualitative inspection of the spectra revealed no ev- 
idence of radial velocity variations associated with bi- 
nari ty, f urther supported by the measurements in Sec- 
tion 



2.2 



Thus, the spectra from each exposure/setting 
were co-added, yielding a signal-to-noise ratio in excess 
of 250 per rebinned pixel from the pipeline. However, 
note that the observations only spanned three nights, 
so our detection sensitivity to radial velocity variations 
from anything but a short-period companion will have 
been low. 

There are archival, low-dispersion observations of 
Skl83 from the International Ultraviolet Explorer (IUE) 
from Hutchings et al. (1991). Although at low resolu- 
tion, the lack of strong UV emission-lines is evident (see 
Figure 4 from Hutchings et al.), consistent with the weak 
wind-profiles seen in other massive stars in the SMC, e.g. 
Walborn et al. (1995, 2000). 

2.1. Spectral Classification 

The AA3950-4750A region of the spectra of Skl83 
is shown in Figure [2 in which the data have been 
smoothed and rebinned to an equivalent resolving power 
of R= 4000 for classification (cf. the new Galactic atlas 
of standards from Sana et al. in preparation). Although 
relatively modest in strength, N IV A4058 emission has 
an intensity comparable to that from N III AA4634-41-42. 
Employing the criteria from Walborn et al. (2002), this 
suggests a classification of 03 V((f*))z, in which the 'z' 
suffix indicates that the A4686 absorption has a greater 
depth than any of the other helium lines (Walborn & 
Blades, 1997; Walborn 2009). 

An early-type classification is supported by the weak 
absorption from Nv (AA4604, 4620), but there are other 
features which are not fully consistent with the 03 clas- 
sification. For instance, the apparent intensities of the 



Setting 


Date 


Exp. Time 


A-coverage 


R 






M 


[A] 




LR02 


2010-10-24 


7x1800 


3960-4564 


7000 




2010-10-26 


2x1800 






LR03 


2010-10-25 


4x1800 


4499-5071 


8500 


HR15N 


2010-10-25 


2x1800 


6442-6817 


16000 



Note. — The final column gives the effective resolving 
power (R) at the central wavelength of each setting. 

He I AA4388, 4471, 4713 lines suggest a later spectral 
type. The intensity of the He II A4026 absorption is also 
relatively strong for an 03 type dwarf, but may be ex- 
plained by a contribution from He I. We note that Hutch- 
ings et al. (1991) classified their low-dispersion optical 
spectrum of Skl83 as 06, but the lack of strong He I 
A4471 absorption in their data argues for an earlier type. 

As a possible explanation of the spectral properties 
of Sk 183 we propose that its spectrum is composite in 
nature, and results either from a binary system or an 
unresolved multiple. While the star appears as a single 
point source in the image from the HST ACS, the spa- 
tial resolution in the Wide Field Channel (2 pixels = O'.'l) 
subtends several thousand AU at the distance of the 
SMC. The composite character of the Sk 183 spectrum 
is further supported by the morphology of the He I ab- 
sorption, which indicates rapid rotation in a secondary 
component (e.g. A4388 in Figure [2]) or perhaps evidence 
for further multiplicity. Without additional data (e.g., 
time-series photometry or multi-epoch spectroscopy) it 
is hard to constrain the nature of Sk 183 further and 
we adopt a classification of 03 V((f*))z + OB. Similar 
composite classifications were employed for three objects 
by Walborn et al. (2002), who demonstrated that one 
(Cyg OB2-22) was indeed a composite of an 03 super- 
giant and an 06 dwarf. 

Circumstantial evidence for the composite nature of 
Skl83 is also provided by the spectra of another star 
in the SMC, MPG 324 (Massey et al. 1989fJ previously 
classified as an 04 dwarf (Niemela et al. 198o; Massey et 
al. 1989; Walborn et al. 1995, 2000; Evans et al. 2006). 
Returning to the FLAMES spectra of MPG 324 from 
Evans et al. (2006), similar features are seen as those in 
Skl83, i.e., N IV emission indicative of an early (<04) 
type, combined with He I absorption suggesting a later 
type (see Figure |2j). Indeed, Evans et al. (2006) noted 
MPG 324 as a single-lined binary due to a ~30kms _1 
shift in the He II A4542 line between two epochs of ob- 
servations with one of the FLAMES settings. The He I 
A4471 line in these two epochs was contaminated by neb- 
ular emission, but the wings of the profile were (qualita- 
tively) consistent with a zero velocity shift, suggesting a 
potentially different origin to the He II absorption. We 
reclassify MPG 324 as 03 V((f*))z + OB, and note that 
Massey et al. (2009) classified MPG 324 as 05 V in 
the course of their atmospheric analysis - they did not 
see N III emission in their (lower-resolution) spectrum, 
but N IV emission remains at A4058, which argues for an 
earlier type. 

9 Aliases: NGC346 #6 (Walborn et al. 1995); NGC346 007 
(Evans et al. 2006). 
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The use of an 03 type to describe both Sk 183 and 
MPG 324 arises from the high-quality of the spectra, in 
which the N IV emission is clearly seen. These classifi- 
cations are based on the spectroscopic framework in this 
regime of massive stars, which have unambiguous crite- 
ria linked to the relative intensities of the nitrogen lines 
(Walborn et al. 2002, 2004; Sota et al. 2011). 

2.2. Stellar Radial Velocities 

Stellar radial velocities (v T ) for Skl83 were deter- 
mined from Gaussian fits to the He II AA4200, 4542 
lines. From measurements of each LR02 observa- 
tion we found weighted means of v r = 162.5 ± 1.2 and 
162.4 ±0.8 km s" 1 from He II A4200 and A4542, respec- 
tively. The quoted errors are 1-er uncertainties on the 
positions of the centroids of the fits. A statistical vari- 
ability analysis of the individual radial velocity measure- 
ments does not detect significant differences, suggest- 
ing a constant velocity from our spectroscopy. From 
fits to the superimposed Ha and [O III] AA4959, 5007 
nebular emission, the mean gas velocity is slightly off- 
set with t; gas = 172.1 ± 1.1 kms -1 . This is approximately 
5kms~ 1 lower than the gas velocity measured at the 
same right ascension (but ~20" further south) by Nigra 
et al. (2008). 

3. PHYSICAL PROPERTIES OF SK 183 
3.1. Distance and Extinction 

The distance to Sk 183 is the largest uncertainty on 
its total luminosity. Given the long-standing questions 
regarding the depth and structure of the SMC, the 'stan- 
dard' distance modulus of ~18.9 mag (e.g. Harries et al. 
2003) might not be appropriate for the Wing. There is 
mounting evidence for objects in the eastern part of the 
SMC to suggest shorter distances, e.g., Howarth (1982), 
Mathewson et al. (1986), Crowl et al. (2001), and Glatt 
et al. (2008). Indeed, in their analysis of the HST imag- 
ing, Cignoni et al. (2009) advocate a distance modulus 
to the young population in NGC602 of 18.7 mag. 

Sk 183 is relatively isolated from bright companions so, 
to determine the line-of-sight reddening, optical UBVR 
photometry was taken from Massey (2002). We also 
employed the F8141F magnitude from Schmalzl et al. 
(2008) and near-infrared JHK S magnitudes from the 
Two Micron All Sky Survey (2MASS; Skrutskie et al. 
2006), as shown in Figure [3] and summarised in Table [2j 
Photometry of Sk 183 at longer wavelengths has become 
available recently from the Wide-field Infrared Survey 
Explorer (WISE, Wright et al. 2010); these data are 
discussed separately in Section [4~2| 

Fr om fits of the adopted model atmosphere (Sec- 
tion 3.2) to the spectral energy distribution of the 
IUE spectroscopy and multi-band photometry (Figure|3| 
we derived E(B - V) = 0.09 ± 0.01 mag. This included 
contributions from both internal SMC and foreground 
Galactic reddening, taking the lower bound of E(B — 
V) = 0.04 mag from Bessell (1991) for the latter. In cal- 
culating the total visual extinction we adopted ratios of 
total- to-selective extinction of 3.1 (Galactic; Howarth, 
1983) and 2.7 (SMC; Bouchet et al. 1985), lead- 
ing to Ay — 0.26 mag and an absolute magnitude of 
My = —5.14 or — 5.34mag (for distance moduli of 18.7 
and 18.9 mag, respectively). Such values are at the faint 



Table 2 

Photometry of Sk 183 



Band 


Magnitude 


rtei. 


u 


12.51 ±0.01+ 


Massey (2002) 


B 


13.59 ±0.01+ 


Massey (2002) 


V 


13.82 ±0.01 


Massey (2002) 


R 


13.89 ±0.02+ 


Massey (2002) 


F814W 


14.128 ±0.003 


Schmalzl et al. (2008) 


J 


14.426 ±0.029 


2MASS 


H 


14.606 ±0.063 


2MASS 


K s 


14.618 ±0.091 


2MASS 


Wl (3.4 fim) 


14.480 ±0.029 


WISE 


W2 (4.6 /an) 


14.255 ±0.041 


WISE 


W3 (12 urn) 


10.368 ±0.043 


WISE 


WA (22 /tm) 


5.955 ±0.031 


WISE 



Note. — t Calculated in quadrature from the errors 
quoted on the colours by Massey (2002). 

end of those found for early-type dwarfs (e.g. Walborn et 
al. 2002), and are comparable to others classified as Vz 
(see, e.g., Figure 4 from Walborn 2009), notwithstanding 
a potential contribution from a companion. 

3.2. Atmospheric Analysis 
3.2.1. PoWR Models 

Physical properties for Sk 183 were estimated from 
comparisons of the FLAMES data with synthetic spec- 
tra calculated with the PoWR code (e.g. Hamann & 
Grafener, 2003; 2004). PoWR treats the wind and the 
photosphere self-consistently and has been used to anal- 
yse the wind characteristics of O-type stars (Oskinova et 
al. 2006, 2007), helium-rich subdwarf O-type stars (Jef- 
fery & Hamann, 2010), and B-type main-sequence stars 
(Oskinova et al. 201 lQ 

In brief, the PoWR code solves the radiative transfer 
equation for a spherically-expanding atmosphere and the 
statistical equilibrium equations simultaneously, while 
also accounting for energy conservation and allowing de- 
viations from local thermodynamic equilibrium (i.e. non- 
LTE) . A PoWR model is defined by its effective tempera- 
ture (T e ff), effective surface gravity (logoff), luminosity 
(L), mass- loss rate (M), terminal velocity of the wind 
(woo), and chemical composition. A standard j3- velocity 
law (/3 = 1) for the wind acceleration was assumed in all 
calculations; the density follows from the velocity via the 
mass-continuity equation. The wind domain is smoothly 
connected to the photosphere, at which point the hydro- 
static equilibrium is approached asymptotically. 

The models included transitions of hydrogen, helium, 
carbon, oxygen, nitrogen, and silicon. Iron and iron- 
group elements were included using the concept of su- 
perlevels, as described by Grafener et al. (2002); this is 
important because of the blanketing effect of the summed 
opacities on the atmospheric structure. 

3.2.2. 'Single-star' Analysis 

We first assumed that the observed spectrum origi- 
nates from a single star. Stellar parameters were de- 

10 Analogous to, e.g., the differences discussed by Rivero 
Gonzalez et al. (2011) between results from CMFGEN (Hillier & 
Miller, 1998) and FASTWIND (Puis et al. 2005), we note there may 
be small systematic differences if the parameters determined from 
PoWR were compared to those obtained with other codes. 
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z K 53 Sk 183 




MPG 324 03 V((f*))z ,+ OB 




Wavelength (A) 



Figure 2. Combined FLAMES spectrum of Sk 183 together with that of NGC 346 MPG 324 from Evans ct al. (2006); both spectra have 
been smoothed and rebinned to an effective resolving power of 4000. Emission lines identified in Skl83 are N IV A4058; Si IV AA4089, 
4116; N III AA4634-41-42, with absorption from N V AA4604, 4620. Absorption lines identified in MPG 324 are He I AA4388, 4471, 4713; 
He II AA4026, 4200, 4542, 4686. Both spectra show a combination of N IV emission and He I absorption inconsistent with a unique spectral 
type, suggestive of composite spectra. The gaps in the spectrum of MPG 324 are due to 'hot' pixels in the detector. 



termined iteratively, with the bolometric luminosity 
calculated fr om the absolute magnitude (—5.14 mag) 
from Section |3.1| The model spectra were convolved 
with a rotational-broadening function to take into ac- 
count an apparent projected rotational velocity of 
vs'mi = 60 ± 20kms _1 , determined from the metal lines 
(N III, N IV, N V, and Si IV). The primary temperature 
diagnostic available is the nitrogen ionization balance, 
with additional constraints given by the Si IV emission 
and the He I A4471 absorption. 

From comparisons of model spectra with the ob- 
servations, we adopted a nitrogen-rich model with 
T e ff = 46±2kK and logoff = 4.0 ±0.1, as shown in Fig- 
ure [3] The stellar radius, at a Rosseland optical depth 
of 20, is R*=9.7R Q . The He II AA4200, 4542 lines were 
somewhat over-predicted by this model. While a cooler 
model (T c ff = 42 kK) still provided a reasonable fit to the 
Si IV emission, it led to significant over-prediction of the 
He I absorption with little improvement for the He II 
lines. A simultaneous fit to the (notably broader) weaker 
He I lines was not possible with a model of a single star. 

The adopted 46 kK model fits the observed N III, N IV, 
and N V features simultaneously with an abundance of 
12+log([N/H]) = 7.73 (±20%). Compared with the Solar 
value of 7.83 (Asplund et al. 2009) this suggests signif- 
icant enrichment compared to the baseline SMC value 



found from H II regions (~6.55, Russell & Dopita, 1990). 
The uncertainty quoted on the T c g estimate (±2 kK) re- 
flects the fitting of the stellar parameters. With more 
modest nitrogen abundances, e.g., a factor of two lower 
(comparable to that obtained for NGC 346 MPG 355 by 
Walborn et al. 2004), the N IV and N V features could 
be fit with a hotter temperature (of a few kK), but a 
simultaneous fit of the N III lines was not possible. 

The mass-loss rate appears to be low, as might be ex- 
pected given the reduced metallicity of the SMC (e.g. 
Bouret et al. 2003). We calculated a grid of models 
where all parameters were fixed, but with a range of 
10~ 9 < M < 10 -5 Mq yr^ 1 . From comparison of the mod- 
els with the observed spectra, M =10 M^yr^ 1 was 
the best compromise to fit the observed UV (N v) and 
optical (Ha and He II A4686) profiles. For example, if 
M > 10~ 7 M Q yr _1 the predicted He II A4686 absorption 
was weaker than observed, suggesting this as an upper 
limit. A more precise estimate of M would require high- 
quality UV and/or if -band spectroscopy. 

To estimate Woo we employed the empirical relation 
from Prinja (1994) between Voo and the wavelength dif- 
ference (AA) between the observed peak emission and 
minimum absorption of the N V or C IV lines in low- 
resolution IUE spectra. Using Prinja's equation (1) 
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for Skl83, AA(N v)«12A, giving Voo w 2150 km s -1 
(± ~300kms -1 ); this value was adopted in our models. 

The degree of clumping in the wind was poorly con- 
strained by the available data because of the low mass- 
loss rate. The final model was therefore calculated with- 
out inclusion for the effects of clumping, but note that 
the exact values of and the clumping factor did not 
have a significant impact on the other physical parame- 
ters. In particular, the N IV emission in our PoWR model 
does not arise from the wind, so was not affected by the 
adopted clumping factor. 

When correcting the effective gravity for the contri- 
bution of radiation pressure we found log g = 4.11, yield- 
ing a spectroscopic mass estimate of M spec = 46lg M . 
This is in good agreement with the evolutionary mass 
obtained from comparisons with the theoretical evolu- 
tionary tracks for the SMC from Brott et al. (2011). 
For usini (init.) < 350 km s , the T e s and L determined 
for Skl83 lies roughly equidistant between the 40 and 
50 Mq tracks, suggesting an evolutionary mass for Sk 183 
of 45 ± 5M , with a corresponding age of ~2 Myr. In the 
LMC grid from Brott et al. there are also models for an 
initial mass of 45 M Q . Scaling the moderate rotation 
models (wsini (init.) < 350 km s -1 ) to the hotter temper- 
atures obtained for stars in the SMC, we find reason- 
able agreement with the parameters obtained for Sk 183 
(Dr I. Brott, private communication). 

The parameters determined from the single-star anal- 
ysis of Skl83 are summarised in Table |3j For complete- 
ness, in the third column we provide the equivalent re- 
sults assuming the canonical distance modulus to the 
SMC of 18.9 mag. 

3.2.3. Composite Model 

We also investigated a scenario in which the observed 
spectrum of Sk 183 is comprised of a hot O-type primary, 
combined with a cooler, lower luminosity, secondary com- 
ponent. Composite spectra were synthesized by scaling 
the continuum fluxes of the two components by their rel- 
ative luminosities, such that the combined model repro- 
duces the observed spectral energy distribution of Sk 183 
(assuming a distance modulus of 18.7 mag). Given the 
large luminosity of the 03 component, the total contri- 
bution from, e.g., a B-type companion, is relatively small 
(indeed, it is within the uncertainty of the single-star fit). 

Adopting the same radial velocity for the secondary 
spectrum, a reasonable fit was obtained with a composite 
of a slightly hotter, early O-type model (T c g = 47.5kK) 
with that of an early B-type giant (with logL/L© =4.1, 
T e fj = 21kK, and log <? e ff = 3.5), as shown in Fig- 
ure [4j A large (projected) rotational velocity of 
v s'mi = 250 km s -1 was necessary for the secondary com- 
ponent to reproduce the broad He I lines. 

While there is only tentative evidence of a secondary 
companion, we include this analysis to explore its poten- 
tial consequences on the combined spectrum. Indeed, if 
a genuine binary, it would be surprising (in an evolution- 
ary sense) to find a B-type giant with a much younger 
03-type dwarf. However, this is not a unique solution for 
the secondary, as hotter models for the secondary (e.g. a 
late O-type star) also give reasonable fits. 

The putative composite model has some advantages 
compared to the single-star fit. In addition to the fits 



to the broad He I lines, the reproduction of the He II 
AA4200, 4542 lines is improved due to dilution of the 
continuum level by the second component. Moreover, a 
more consistent fit is achieved for the He I A4471 line - 
we were unable to reproduce the broad wings and nar- 
row core of this line simultaneously in the single-star 
model (see Figure [5]). The parameters of the primary 
and secondary components in the composite model are 
summarised in Table [3] 

4. DISCUSSION 

The discovery of such an early- type star in the SMC, 
let alone in the Wing, is noteworthy. To place its rar- 
ity in context, there are published classifications for over 
5000 luminous stars in the SMC (e.g. see compendium 
by Bonanos et al. 2010) but only two with secure spec- 
tral types of earlier than 04, namely NGC 346 MPG 355 
and AzV 435^) Thus, including the revised classification 
of MPG 3247~the new spectroscopy of Sk 183 brings the 
total of 02 and 03 stars known in the SMC to only four 
objects. 

These stars provide a glimpse of some of the most 
massive O-type stars in a very metal-poor environment, 
into a domain of significantly diminished mass-loss rates 
(e.g. Bouret et al. 2003). Indeed, in their discussions 
of chemically-homogeneous evolution in the context of 
gamma-ray bursts, Yoon et al. (2006) noted NGC 346 
as the only obvious cluster in which to examine the evo- 
lution of the most massive stars in a truly metal-poor 
regime; Sk 183 provides an excellent additional case for 
further study. In particular, spectroscopic monitoring 
would be helpful to unveil its true nature. 

4.1. Ionisation of the N90 Nebula 

The Spitzer observations revealed a population of 
young stellar objects along the inner edge of the N90 neb- 
ula, with Carlson et al. (2007) concluding that star for- 
mation started in the central cluster ~4 Myr ago and has 
propagated outwards since then. The peak star forma- 
tion rate in the region has occurred over the past 2.5 Myr 
(Cignoni et al. 2009), during which time it has been ap- 
proximately constant (Carlson et al. 2011). The evolu- 
tionary age obtained for Skl83 (~2Myr, Section 3.2.2), 
is consistent with it forming in this recent phase of star 
formation. Meanwhile, the gas in the interior of the H II 
region is nearly quiescent and does not show violent mo- 
tions (Nigra et al. 2008), leading Carlson et al. (2011) 
to conclude that photoionisation is the primary trigger 
of the current star formation. 

N90 has been observed from the ground and with 
the HST in narrow-band filters, thus isolating the Ha 
and [N 11] emission. Using the ACS images obtained 
in program GO-10284 (P.I.: A. Nota) we integrated 
the flux from the high surface-brightness component of 
the nebula finding ^2 x 10~ n ergs -1 cm -2 , after correc- 
tions for the foreground extinction and the weak [N 11] 
emission. The former was taken as Ar = 0.16 mag from 
the NASA/IPAC Infrared Science Archive Extinction 

11 MPG 355 was classified as 02 III(f*) by Walborn et al. (2002) 
and ON2 III(f*) by Walborn et al. (2004); AzV 435 was classified 
as 03 V((f*)) by Massey et al. (2005). Note that the spectral type 
listed by Bonanos et al. for AzV 14 is 03-4 V (Garmany et al. 
1987), but this was revised to 05 V by Massey et al. (2004). 
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Figure 3. Upper panel: model fit to the spectral energy distribution of Skl83, combining IUE spectroscopy and multi-band photometry 
as indicated. The total linc-of-sight reddening, taking into account Galactic and SMC contributions, is E(B — V) = 0.09 mag. Middle and 
lower panels: comparison of the optical FLAMES spectrum of Skl83 (blue solid line) with the adopted 'single star' PoWR model (red 
dashed line). 



Calculator (which uses the results from Schlegel et al. 

1998) , while a contribution to the [N II] emission of 4% 
was adopted (e.g., Peha-Guerrero et al. 2012). How- 
ever, the ACS image covers only the core part of N90, 
which extends over a larger area of 8' x 6'. Thus, inte- 
grating over the entire extent of N90 from the Magel- 
lanic Cloud Emission Line Survey (MCELS; Smith et al. 

1999) we obtained / (Ha)«8x 10" 11 ergs" 1 cm" 2 (cf. 
/ (Ha) = 10.8 ±1.1 x 10" 11 ergs" 1 cm" 2 from Kennicutt 
& Hodge, 1986). The flux from the MCELS image corre- 



sponds to a total luminosity of > 3.2 ± 0.3 x 10 37 ergs -1 
(with the uncertainty indicative of the two distance mod- 
uli, i.e. 18.7 and 18.9 mag). 

Converting this L(Ha) to a Lyman continuum lumi- 
nosity for a standard nebular case yields an ionization 
rate of Qo « 2 xlO 49 s . This will be a lower limit 
as some ionization surfaces within N90 are hidden from 
optical view by dark clouds and we have not accounted 
for leakage. Aside from Skl83, the next earliest objects 
known in NGC 602a are stars #5 and #2 (Hutchings et 
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Table 3 

Derived Properties of Sk 183 





Single-star Model 


Composite 


Model 




DM = 18.7 mag DM = 18.9 mag 


Pri. 


Sec. 




—5.14 mag —5.34 mag 


— 5.14 mag 


log L/L/7S 


5.58 ±0.02 5.66 ±0.02 


5.57 


4.1 




46 ± 2 kK 


47.5 kK 


21 kK 


R* 


9.7-Rq 10.7R© 


9.6 Rq 


8.5 Rq 


log g e ff 


4.0±0.1 4.0±0.1 


4.0 


3.5 


log g 


4 11-0.07 4 n -0.07 


4.11 


3.5 


-^spcc 


46± 9 S Mq 541 11 M 


40t* Mq 


8.7 Mq 


M 


10~ 7 Mq yr- 1 


10" 7 Mq yr" 1 




v sini 


60 ±20 kms -1 


60 km s -1 


250 kms" 1 


Voo 


2150±300kms- 1 


2150 kms- 1 




log Qo 


49.37s- 1 49.45s- 1 


49.41s- 1 




12+log([N/H]) 


7.73 (±20%) 


7.73 





Note. — The quoted errors on the single-star models are from fitting uncertain- 
ties. We do not quote uncertainties on the composite model as it is non-unique fit 
to the data. 



al. 1991, classified as 08 and 09, respectively). The ion- 
ising fluxes from late O-type stars are over an order of 
magnitude lower (see, e.g., Smith et al. 2002), such that 
Sk 183 appears to be the dominant source of ionisation. 

4.2. Mid-infrared Excess 

The spectral energy distribution of Sk 183 is well 
matched in the near-IR by both the single-star and com- 
posite models (Figures [3] and PO) . However, inspection of 
the new WISE photometry (Table [2]) reveals an inter- 
esting upturn in the flux distribution, as shown in Fig- 
ure [m Observations with the short-wavelength WISE 
bands (Wl and W2, corresponding to 3.6 and 4.6 /im, 
respectively) are in good agreement with the predicted 
flux from our PoWR models. But at the longer wave- 
lengths of the W3 and WA bands (12 and 22 /im), there 
is a clear mid-IR excess compared to the predicted distri- 
bution. A peak of the mid-IR flux distribution at 22 /im 
would argue for a blackbody temperature of order 130 K, 
or even cooler if the flux peaks at longer wavelengths. 

Similar mid-IR excesses for otherwise normal OB-type 
stars have been seen in the SMC (Bolatto et al. 2007; 
Ita et al. 2010; Bonanos et al. 2010) and LMC (Evans 
et al. 2011). To investigate the mid-IR excess in Skl83 
further, we obtained images of the region from the Spitzer 
SAGE-SMC Survey (Gordon et al. 2011). 

The InfraRed Array Camera (IRAC, 3.6, 4.5, 5.8, 
8 /im) and Multiband Imaging Photometer for Spitzer 
(MIPS, 24 /xm) images of Skl83 are shown in Figure [7l 
together with an Ha image obtained with the M7> 
saic II camera on the 4-m Blanco telescope at the Cerro 
Tololo Inter-American Observatory. From photometric 
extractions of the IRAC images we found: [3.6] = 14.75, 
[4.5] = 14.71, [5.8] = 14.67, and [8.0] = 14.66 mag. How- 
ever, the peak of the 24 /im flux is offset slightly and 
appears extended (cf. the Airy rings seen around the 
point sources to the east), suggesting it as nearby dust 
emission, or unresolved dusty shocks (e.g. those seen in 
the Carina nebula, Smith et al. 2010), but not directly 
associated with Sk 183. 

4.3. High-mass Star Formation in NGC 602 

In the standard model of star formation, the initial 
mass function (IMF) is purely statistical in nature and 



the probability of very massive stars being formed, while 
low, is independent of the mass of the star-forming re- 
gion (e.g. Parker & Goodwin, 2007). In the integrated 
galactic IMF model (IGIMF, e.g., Weidner & Kroupa 
2006) this probability increases with the mass of the 
star-forming region, with the prediction that the most 
massive stars should be found in the largest star-forming 
complexes. 

NGC 602a is a relatively low mass cluster, with a stel- 
lar mass of ^2000 M Q derived from the combination of 
the optical (Cignoni et al. 2009) and infrared imag- 
ing (Carlson et al. 2011). In this context, the IGIMF 
predictions from Weidner et al. (2011) suggest a maxi- 
mum stellar mass of ~50 to 60 Mq, roughly consistent 
with the estimated mass of the primary of Sk 183. In- 
triguingly though, Sk 183 is located approximately 45" 
to the north-west of NGC 602a (corresponding to a pro- 
jected distance of 12-13 pc, depending on the adopted 
distance modulus). This suggests the star has formed 
in relative isolation compared to the main body of the 
cluster (with no obviously massive cluster visible within 
3pc, e.g., Bressert et al. 2012), or that it is a slow run- 
away from NGC 602a, as alluded to by Gouliermis et al. 
(2012). Analysis of the line-of-sight velocities of other 
stars in N90 (and the gas profiles super-imposed on their 
spectra) is underway and should provide a more complete 
picture of the dynamics of this region. 

While the single case of NGC 602 can not distinguish 
between the two IMF models, the SMC is rich in small, 
star-forming regions with masses in the range of 1000- 
2000 Mq. A more complete census of O-type stars in 
SMC clusters would provide an interesting test of the the- 
ories of the high-mass IMF, complementing recent work 
investigating 'isolated' O-type stars in the SMC by Lamb 
et al. (2010) and Selier et al. (2011). 
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Figure 7. Mosaic II (Ho), Spitzer IRAC, and MIPS (2'x2') images of Skl83 (marked by the ope n cro ss) and the surrounding nebula. In 
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